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We study the fate of a localized wavepacket in a classical conformal field theory with attractive 
interaction V(<t>) — — \ A(/> 4 . As potential is unbounded from below, homogeneous field collapses to 
singularity in finite time. However, finite size wavepacket can disperse before it collapses. Competi- 
tion between the two outcomes results in a critical behavior, much like the one seen in gravitational 
collapse. We calculate the critical exponents, and show that there are static regular soliton-like 
solutions in the theory. 
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I. INTRODUCTION 

Scale-invariance is ubiquitous in physics, appearing 
in statistical mechanics and phase transitions, quan- 
tum field theories near fixed point, and the much-lauded 
AdS/CFT correspondence [T]. In this paper, we explore 
dynamics of a simple classical conformal field theory in 
a fiat spacetime, consisting of a single real- valued scalar 
field with action 
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and the quartic potential with the "wrong" sign 
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corresponding to attractive self-interaction. Although 
potential is unbounded from below and the theory is 
unstable, the classical vacuum <j> = is infinitely long- 
lived, and the corresponding quantum theory is renor- 
malizable [2j [3]. In fact, this model is in some sense 
better off than familiar gravity, which is also attractive, 
also has unbounded negative binding energy, singulari- 
ties appearing in evolution of a regular initial data, but 
is non-renormalizable on top. 

Quantum field theory with qb A interaction has been ex- 
tensively studied, with renormalization group calculated 
perturbatively to fifth loop long time ago [41. Recently, 
Dvali et. al. [5] argued that the theory (l]2) captures 
some features of quantum phenomena, asymptotic free- 
dom and dimensional transmutation, at a classical level. 
This would be pretty interesting, because in classical field 
theory we have direct access to fully non-perturbative 
and non-linear solutions, which are readily found using 
numerical methods. Models of the type Q might also 
have cosmological applications, where conformal symme- 
try can be used to generate scale- invariant primordial 
fluctuations [6j[7]. 



With this in mind, it would seem that the theory ([2| 
warrants some further investigation. Seeming simplicity 
of the model is deceptive, and there are some rather non- 
trivial things about it, which is what I report here. 

Evolution of the classical scalar field with potential 
Q is governed by the wave equation 



□0 + A0 3 = 0. 



(3) 



Conformal symmetry can be used to make variables di- 
mensionless by rescaling to an arbitrary length scale t. 
In addition coupling A can be rescaled to unity by 



(4) 



* Electronic address: frolov@sfu.ca 



as long as it does not vanish, which is what we will do in 
most calculations below. 

The outline of the paper is as follows: In Section [ll| I 
explore the eventual fate of a localized wavepacket and 
find that there is dynamical critical behavior in this the- 
ory, similar to the critical gravitational collapse [8]. In 
Section III, I derive static spherically symmetric solu- 
tions, which turn out to have regular localized soliton- 
like states that are unstable and probably correspond to 
critical solutions of the dynamical collapse. As the en- 
ergy of these is infinite due to slowly decaying tails, I 
construct static non-linear dipole configurations in Sec- 
tion [TVl and show that their energy is finite. I conclude 
with some speculations in Section [V] Finer points of clas- 
sical beta function asymptotics and behavior of irregular 
static solutions are hidden in Appendix [Al 



II. CRITICAL COLLAPSE 

Conformally invariant theory with negative potential 
will have a direction in which the potential is unbounded 
from below, and hence is unstable. This is not necessar- 
ily fatal by itself. As long as the vacuum does not decay 
explosively, this instability simply means that fields are 
prone to collapse under self-attraction, just like Jeans 
instability in gravitating matter. Just as in gravity, sin- 
gularities can form as a result of a field collapse in an 
attractive CFT. 
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FIG. 1: Dispersal of a momentarily stationary spherical Gaussian wavepacket in free theory (a), and in no n- linear CFT with 
initial amplitude slightly below (b) and slightly above (c) critical value. In the latter case, singularity forms in finite time. 
Density plot in {r, £}-plane shows isosurfaces of q(r, t) = r0, with evolution of field profiles 0(r, t = const) overlaid in pale 
shading. Solid red line plots field value at origin 0(r = 0,£) as a function of time. 
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FIG. 2: Scattering of an ingoing sine squared wavepacket in free theory (a), and in non-linear CFT with initial amplitude 
slightly below (b) and slightly above (c) critical value. In the latter case, singularity forms in finite time. Same legend as Fig.[l] 



So, what is the eventual fate of a given field distribu- 



tion in an attractive CFT? 
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(a) (b) 
FIG. 3: Scaling of lifetime r of momentarily stationary (a) and scattered (b) wavepacket with initial field amplitude p. 



Dynamics of the homogeneous field evolution are the 
same as for a particle falling down in potential ([2]); a 
particular solution of Q is 



<f>(t) 



t/U ' 



2 1 

A 0o* 



(5) 



The field collapses from finite value </>o to singularity in 
a finite time £*, with the timescale set by the field value. 
The vacuum solution, although unstable, takes infinitely 
long to decay as the potential is rather flat at the top. 

Localized wavepacket is not going to hang around wait- 
ing for the collapse to happen. Characteristics of the 
wave equation ([3| propagate at the speed of light, so the 
initially stationary wavepacket will disperse in a time 
roughly comparable to its size. Competition between 
propagation and self-attraction results in a critical be- 
havior much like the one seen in gravitational collapse 
[8HTU], but in a much simpler theory. 

To analyze what happens, lets solve the equation (|3| 
for a spherical wave 



dt 2 



d 2 (j) 
dr 2 



2 d<j> 
r dr 



Xcjf = 0. 
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With the usual variable change 0(r, t) = q(r,t)/r, it re- 
duces to a one dimensional wave equation 



dt 2 



d 2 q x q 3 
+ A 







with Dirichlet boundary condition at the origin r 







r= 



dq 
dr 



dq 
dt 
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supplemented by outgoing boundary condition at outer 
domain boundary r = r c for simulation. The semi-linear 
partial differential equation is straightforward to solve 
numerically; we use second order leapfrog method on a 
large uniform grid (2 15 nodes for results quoted here). 



Evolution of initially stationary Gaussian wavepacket 
4>(r,t = 0) = pexp(— r 2 / s 2 ) is shown in Fig. [l] In free 
theory ([l^), momentarily stationary configuration splits 
into a left and right moving waves, both of which dis- 
perse to infinity after the left mover reflects off of origin. 
For weak non-linearity, the evolution is pretty much the 
same, but as one increases the amplitude p (keeping s 
constant), situation changes. The wavefront comes out 
without the phase reversal, leaving a blob of scalar field 
behind, which is held by attractive interaction. If the am- 
plitude is slightly smaller than a certain value p*, it hangs 
close to balance but disperses eventually as in Fig. (TJd. If 
one increases the amplitude even slightly above critical 
value p*, singularity forms as in Fig. ^p. The change is 
very abrupt, the difference of initial amplitudes between 
Figs.[ljD andP is roughly 0.1%. Same thing happens in 
a scattering of initially ingoing wave, as shown in Fig. [2j 

This is a by- now- familiar picture of the critical gravi- 
tational collapse [8HTU] , except it happens in a pure field 
theory in a flat spacetime. One can take inverse lifetime 
1/r (measured from the moment the maximum of the 
pulse arrives to the origin to the moment singularity first 
forms) as an order parameter for the phase transition. It 
is zero in a regular phase, and scales with amplitude as 
shown in Fig. [3] in a singular phase. The scaling is a bro- 
ken power law: For large amplitudes scaling is linear as 
characteristic of homogeneous collapse ([5|, but close to 
criticality it switches over to much smaller exponent with 
a numerical value of approximately 1/16, which is suffi- 
ciently close for different families of solutions to suspect 
some sort of universality. If the scaling persists as one 
tunes the parameter to criticality, the critical solution 
separating the two phases would be infinitely long-lived, 
i.e. quasi-static. This is unachievable in collapse simula- 
tion, but is easy to find directly if it exists, which we will 
do next. 
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III. STATIC SOLUTIONS 



As we have seen in the last section, arbitrarily long- 
lived configurations can be obtained in the field collapse 
by fine-tuning a single parameter. In fact, the theory 
admits one-parameter family of regular soliton-like solu- 
tions. Spherically symmetric static solutions satisfy 



dr 2 



2 dcj) 
r dr 



\<tf = 0. 



(9) 



Regularity demands that the power series expansion of 
the solution around the origin is an even function of r, 
while scaling symmetry guarantees that variables enter 
the expansion only in a combination £ = Xcj^r 2 



(r) * 4>o 
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(10) 



Unlike free theory, solutions with A > will decay at 
infinity regardless of the value 0o field takes at the origin. 
The easiest way to see this is to do a variable change 



(j>{r) 



(3r)l 



3r 



(ii) 



which casts the equation ([9| into the form of time- 
dependent anharmonic oscillator 



0. 



(12) 



The above equation corresponds to canonical equations 
of motion for an (x-dependent) Hamiltonian 



1 



^ 2 



1 



(13) 



which makes it clear that any solution will go to a limit 
cycle at large x. Scaling A and amplitude of oscillations to 
unity, the limit cycle solution is given in terms of Jacobi 
elliptical cosine function, or its harmonic expansion [11] 

l 3 v~~^ cos2(?7/ — h)kx , 

i>(x)=cn(x,2-*)=2-ikY, ^ T\ ' (I 4 ) 

^-^ cosnm — ^ )7T 

n=l v ^ y 

It is periodic with period 2ir/k = 7r~^r 2 (|), and the 
harmonic expansion is exponentially converging, so only 
a few terms are needed to accurately represent its shape. 
While Hamiltonian (13) is integrable, its integrals of 



motion are not easily obtainable in closed form [I2j [13] , 
so we find the solutions by numerical integration. Two 
initial conditions are required to specify the solution of a 
second order differential equation; one is the field value 
at origin 0o, the second is supplied by regularity condi- 
tion (10). Thus we have a one-parameter family of static 



everywhere regular solutions, which are localized around 
origin as shown in Fig. [4] and oscillate and decay away 



as r 



-2/3 




FIG. 4: Static spherically symmetric solutions (j)(x) in non- 
linear CFT as a function of rescaled radius x = (Sr) 1 ^ 3 . 
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FIG. 5: Regular static spherically s ymm etric solution <j)(x) in 
green versus global approximation (16) in red. 



at spatial infinity. 



A closer inspection of numerical results reveals that 
the regular solution transits from its behavior at origin 
to the limit cycle almost exactly like a cardinal function 

il> c (x) = sn(x, 2 - ^) - cn(x, 2~i). (15) 

Swapping elliptical functions for trigonometric and in- 
troducing a few more coefficients to match power series 
at origin (10), we obtain a global approximation of the 
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FIG. 6: Bipolar coordinate system. Lines p — const (shown 
in red) and — const (shown in green) are orthogonal. 



regular solution in terms of elementary functions only 



sin kx \ > 

t/>(x) ~ ( -cos /cx J ^a n 

n=0 



kx 



2n 



sm(kx) 
(kx) 2n 



(16) 



where coefficients a n are selected to keep only powers of 
x 6 in expansion around origin. As Fig. [5] shows, approx- 
imation ( [T6| ) reproduces the exact solution obtained by 
numerical integration remarkably well. 

Solutions that blow up at origin require point-like ex- 
ternal charges to source them, but existence of solitons 
implies that the value of the source does not determine 
the solution uniquely. To illustrate the point, Fig. [4] also 
shows static solutions with positive and negative diver- 
gence at origin "riding" on top of the soliton and having 
virtually the same asymptotic away from the origin. Be- 
havior of irregular solutions is substantially more com- 
plex then assumed in [5 . As explained in Appendix \K\ 
it involves two distinct length scales, not one. 



IV. DIPOLE BOUND STATES 

Although the solutions described above are everywhere 
regular and localized around the origin, the field tails 
decay too slowly (</> ~ r ~ 2 / 3 a t large r) for the total 
energy of the state to be finite. This is never a problem 
in practice as you always deal with finite truncations if 
you are trying to assemble the state dynamically, but is a 
reason we put a qualifier "soliton-like" instead of calling 
it a soliton. In this section, we show that in the soliton 
- anti-soliton pair the slow-decaying tails interfere away, 
and the total energy is finite. Such a dipole pair forms 
an unstable bound state. 



We will look for the static dipole solutions by spectral 
decomposition in bispherical coordinate system, which is 
adapted to the symmetries of the solution we seek [13] . 
Bispherical coordinates are obtained by axial rotation of 
a bipolar reference frame, which in turn is a conformal 
transformation of a (p, #)-plane by a factor 



Q — cosh p — cos 0, 



(17) 



which leaves coordinate grid orthogonal, but shrinks in- 
finity p = ±oo to two points - the location of dipole foci 
- as illustrated in Fig. [6j Explicit transformation back 
to Cartesian coordinates reads 



x 

y 



ft 1 sinh p, 

sin#cos(/?, 

ir 1 - 



(18) 



suit/ sin (/?, 

while the flat space metric in bispherical coordinates is 

ds 2 = Q~ 2 (dp 2 +d6 2 ±sin 2 6> dip 2 ). (19) 

Bispherical coordinate system allows separation of vari- 
ables for Laplace equation, so we will decompose the so- 
lution we seek onto a countable set of basis functions 

<j)(p, 0,<p) = ^ U k emFk£m(p, if) (20) 
him 

formed by a tensor product of a radial basis Rk(p) and 
the usual spherical harmonics Ye m (0,(p) 

F Mm (p,0,f) = (2Sl)*R k (p)Y ern (9,<p). (21) 

With a little work, one can show that the variables in 
Laplace equation indeed separate, with 

AF Mm = 2^q! [R'l - i(2*+ l) 2 R k ] Y £m . (22) 

The natural choice for radial basis in dipole symmetry is 



Rk(p) = sinh (h + \)p, 
which makes the Laplace operator diagonal 



AF k £m = ^ 



(k 



Fk£n 



(23) 



(24) 



and the free field dipole solution formed by a pair of point 
charges +q and — q simply 



^dipole — qFooo(p,Q)- 



(25) 



However, as exponentials are rapidly growing functions 
as p — » ±oo, the convergence rate for coefficients of ex- 
pansion in this basis is abysmal. Instead, we notice that 
basis functions (23) are simply odd power polynomials 



in sinhp/2, and recombine the basis set into Chebyshev 
polynomials of argument rescaled inside a bounded do- 
main \p\ < 

Pmax 



Rk(p) = T 2 k+i (crsinh , a 1 = sinh - 



(26) 
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FIG. 7: Dipole solution in free theory (a), and a sequence of non-linear static solutions (b,c,d) as amplitude of the field around 
focal points is increased. Shaded contours show various isolines = const, with thick white one highlighting = 0. Empty 
regions around focal points are excised from the solution domain, and might require charge distributions for regular completion. 



The simplest way to evaluate the new radial basis is by 
using trigonometric representation 

Rk=cos(ng), cos g = a sinh ^ , n = 2k + 1, (27) 

while the second derivative of R^ with respect to p which 
enters the Laplace operator can be evaluated as 



R 



,, 1 sin no , 9N 1 cos no 9/ 9 9N 

— _^ n (l + a 2 ) cos £ - - - - 2 2 - 1 
4 sin q 4 sin" £> 



2 / 2 ,2 
n (^cos £> + a 



(28) 

Evaluated at a set of points {a? a }, Laplacian and field 
value operators ( |20|21 ) are just big matrices acting on a 
vector of spectral coefficients u drawn from a linear space 
spanned by values of composite index him in ( 20 ) 



1 • u, (j)(x) = F • u. 



(29) 



Truncating the spectral expansion and selecting an ap- 
propriate collocation grid to evaluate the residuals, linear 
elliptical partial differential equation A0 = f(x) turns 
into a linear algebra problem D • u = / [15 . An optimal 
grid is determined by the basis choice, which we take to 
be the Gauss-Lobatto grids for Chebyshev and Legendre 
polynomials in radial and angular directions, respectively 



7T TV- 



Qi 



2N- 



PN-licOsOj) = 1. 



(30) 



Direct inversion of the resulting linear algebra problem 
solves linear elliptic PDE to arbitrary precision [T5] . 
Alas, the same cannot be said even for semi-linear PDE. 

It is well known in computational mathematics that el- 
liptic semi-linear PDEs like A0 + 3 = admit infinitely 
many distinct solutions even in compact domains, even 
with trivial boundary conditions [16]. Without any ad- 
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ditional constraints, positions and amplitude of the soli- 
tonic waves will be undetermined. To break this degen- 
eracy, we impose explicit Dirichlet conditions on <p at ar- 
bitrary excision boundary p = p max close to focal points, 
and look for the solution outside. As all of the basis func- 
tions decay as 1/r 2 , no additional boundary conditions 
arise at spatial infinity (p = 0, 6 = 0). Then, starting 
with some initial guess, one tries to zero the residual 



b(x) = A^) + A0 3 -/(x) = 



(31) 



by an iterative scheme. As one soon finds out, Newton- 
Kantarovich iteration obtained by linearizing the above 
equation 



(D + 3A0 2 F) • Su = -b 



(32) 



fails to converge if non-linearity is sufficiently strong (and 
it is for oscillatory solutions we seek). So one is re- 
duced to crawling down the direction of steepest descent 
in Levenberg-Marquardt fashion 



(l t l) • Su = -L T • B, 



(33) 



and hoping for the best. Here tilde denotes the usual 
Levenberg-Marquardt regularization operator scaling the 
diagonal elements of a matrix by a factor of 1 + a, and L 
and B are the preconditioned Newton iteration matrices 



L = fri (B + 3A(/> 2 F) 



B = ft"* b. 



(34) 



Preconditioning by Q 2 reflects the volume factors in the 
residual squared cost function 



X 



I 



b 2 (x)Vgd 3 x. 



(35) 



Thus, after much effort and some applied magic having 
to do with initial guess, the solutions shown in Fig. [7] 
were obtained. As degree of non-linearity is increased, 
the solution goes from a simple deformation (b) of a free 
dipole (a) through a sequence of sign reversals (c,d) which 
change the topology of <\> = surface. All of the dipole 
solutions have finite energy, as far away from the dipole 
field (j) decays as 1/r 2 . 



V. DISCUSSION 

The most interesting result of this paper is that a ver- 
sion of critical collapse exists in a theory without gravity, 
and a simple one at that. The criticality is driven by an 
attractive self-interaction of the field. The critical so- 
lutions are long-lived, and are likely the static soliton- 
like solutions we found. Inverse lifetime of a collaps- 
ing wavepacket scales linearly with amplitude for large 
amplitudes or homogeneous configurations, but switches 
over to a much shallower power law close to criticality. 
The value of critical exponent is determined by the spec- 
trum of unstable perturbation of the critical solution [9] . 



This we leave for future study, but since approximate uni- 
versality is observed in the critical collapse of arbitrary 
initial data, the spectrum should have a mass gap, or be 
discrete. 

We further investigate static spherically symmetric so- 
lutions. The theory admits everywhere regular localized 
vacuum solutions - solitons, but their energy is infinite 
due to slow decay of the oscillating tails at spatial in- 
finity. Irregular solutions have an asymptotic at origin 
where effective coupling constant runs logarithmically, 
as the one in quantum chromodynamics. This, taken 
together with infinite energy of an isolated charge, lead 
to an interpretation of this theory as showing a classical 
analogue of confinement [5,. Leaving aside the matter 
of solutions being unstable, we find the story is substan- 
tially more complicated — in a general static solution two 
other intermediate asymptotic regimes interject between 
logarithmic running and the limit cycle contributing to 
energy divergence. 

We also construct static non-linear solutions for soliton 
- anti-soliton pair. Such a dipole has finite energy, and 
thus forms a bound state. As a final note, we must re- 
mark that the classical conformal field theory we consider 
here is deceptively simple, but actual non-linear solutions 
have very rich phenomenology. In many ways they are 
similar to solitonic solutions of Korteweg-de Vries equa- 
tion and Petviashvili monopole vortices encountered in 
atmospheric science [17] . If one breaks conformal invari- 
ance by introducing mass terms, static solitons discussed 
here will probably turn into oscillating breather modes, 
not quite like in [T5J HH] , but rather related to oscillons 
discovered in [20, 2T] . 
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Appendix A: Classical Beta Function 

As the theory is scale- invariant, the equations of motion 
can be rewritten in an autonomous form. Introducing variable 
q = rcj) once again and switching over to logarithmic spatial 
scale j! = lnr, equation (|9l becomes [5] 



2 - q + Xq = 0. 



(Al) 



Dynamical system can be analyzed in terms of classical renor- 
malization group variables, effective coupling and its running 



a(n) = \q 2 (p), P(n) = a(p). 



(A2) 



As with any autonomous dynamical system, the order can be 
reduced by considering momentum a function of generalized 
coordinate /3 — /3(a), not the scale /x 



{? = {,+ f -2c? ={}£. 

2a da 



(A3) 
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FIG. 8: Classical renormalization group flow in terms of 
a = Xq 2 and f3 = da/ din r for static spherically symmetric 
solutions. Thick green line is a regular solution of Fig. [4] red 
and blue lines are irregular solutions with positive and nega- 
tive divergence at origin, black line is asymptotic /3 = 2a 2 . 



Behavior of solutions as r — » cor resp onds to an UV fixed 
point at a = 0, but the equation ( |A3[ ) admits two distinct 
flows to the fixed point, (3 = 2a and f3 = 2a 2 (and only two 
in the power law class) . Substituting ansatz f3 = ca 1 in ( A3 ) 



-(1 - 2j)a 



2 7 -l 



+ ca 1 



2a z = 0, 



(A4) 



one can see that for leading terms to cancel as a — >• 0, the 
power exponent in the second term has to coincide with the 
first or the last one. First alternative yields /3 = 2a, the 
second ft — 2a 2 . Integrating renormalization group equations 
give corresponding asymptotic solutions 



2a 



■ 2a 



a 2 2 

Aq (X r , 



■ 2 In 



const, 



(j) oc 



21n- 



(A5) 
(A6) 



So, which one is chosen by the static solution? 

Turns out the general solution follows both asymptotics in 
turn, abruptly switching over at a certain scale £*, as illus- 
trated in Fig. [8] This can happen because, while the theory 
is scale- invariant, the solutions are not. Regular soliton solu- 
tions are characterized by a single length scale £o = A -1 / 2 ^ 1 . 
Irregular solutions introduce a second length scale £ q = qcf)^ 1 
related to the value of charge sourcing them. Starting from 
r = oo to r — 0, the general solution then goes through four 
disti nct a symptotic regions: the limit cycle (14), the soliton 
core (A5), the transition over q ~ const branch, and finally 



the asymptotically free regime with logarithmic running ( A6 ) . 
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